"Soft" Anharmonic Vortex Glass in Ferromagnetic Superconductors 
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Ferromagnetic order in superconductors can induce a spontaneous vortex (SV) state. For external 
field H = 0, rotational symmetry guarantees a vanishing tilt modulus of the SV solid, leading to 
drastically different behavior than that of a conventional, external-field-induced vortex solid. We 
show that quenched disorder and anharmonic effects lead to elastic moduli that are wavevector- 
dependent out to arbitrarily long length scales, and non-Hookean elasticity. The latter implies that 
for weak external fields H, the magnetic induction scales universally like B{II) ~ B{0) + cll°', with 
a Ri 0.72. For weak disorder, we predict the SV solid is a topologically ordered vortex glass, in the 
"columnar elastic glass" universality class. 
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Rare-earth borocarbide materials exhibit a rich phase 
diagram that includes superconductivity, antiferromag- 
netism, ferromagnetism and spiral magnetic order. Q-]^] 
In particular, there is now ample experimental evidence 
that, at low temperatures, superconductivity and fer- 
romagnetism competitively coexist in ErNi2B2C com- 
pounds. Other possible examples of such ferromagnetic 
superconductors (FS) are the recently discovered high 
temperature superconductor Sr2Y Rui_a;Cua;06 and the 
putative p-wave triplet strontium ruthenate supercon- 
ductor, Sr2Ru O4, which spontaneously breaks time re- 
versal symmetry. For sufficiently strong ferromagnetism, 
such FS's have been predicted 1^ to exhibit a spontaneous 
vortex (SV) state driven by the spontaneous magnetiza- 
tion, rather than by an external magnetic field H. The 
novel phenomenology of the associated SV solid is the 
subject of this Letter. 

Here we will show that for H = 0, the elastic prop- 
erties of the resulting SV solid differ dramatically and 
qualitatively from those of a conventional Abrikosov lat- 
tice. The key underlying difference is the vanishing of 
the tilt modulus, which is guaranteed by the underly- 
ing rotational invariance (but see below). Although this 
invariance is broken by the magnetization, M, the tilt 
modulus remains zero because this breaking is sponta- 
neous. This contrasts strongly with a conventional vortex 
solid, where the rotational symmetry is explicitly broken 
by the applied field H. All of our conclusions, e.g., the 
unusual B{H) relation, illustrated in Fig.Q are a direct 
consequence of this important observation. 

In particular, we find that this "softness" (i.e., vanish- 
ing tilt modulus) of the SV lattice drastically enhances 
the effects of quenched disorder. As in conventional 
vortex lattices [Q, any amount of disorder Ay, how- 
ever weak, is sufficient to destroy translational order in 
SV lattices. Here the finite ordered domains are diver- 
gently anisotropic, with dimensions oc and 
oc These lengths are measurable in scatter- 



ing and transport measurements. Unlike conventional 
lattices, however, in SV lattices the disorder also qualita- 
tively alters the elastic behavior at long distances, leading 
to "anomalous elasticity" : a universal scaling of elastic 
moduli with wavevector k out to arbitrarily long length 
scales, with some elastic moduli vanishing, and others di- 
verging, as wavevector k 0. Q This behavior is charac- 
teristic of a new kind of topologically ordered "columnar 
elastic glass" (CEG) phase of vortices, which is stable, for 
weak disorder, against proliferation of dislocations. 
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FIG. 1. The nonlinear and universal power-law B{H) 
scaling, which at weakest fields H < Her. and strongest 
fields H > Hnl is cutoff by the crystal symmetry breaking 
anisotropy and (,nl, respectively. 

The best way to experimentally probe a spontaneously 
broken symmetry is to break it directly with an external 
field. We predict that as a consequence of the anoma- 
lous elasticity, the increase in the magnetic induction, 
6B{H) = B{H) — B{0) over the .spontaneous induction 
B{0), due to a weak applied field H along M obeys a 
universal "non-Hookean" scaling law: 
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with the universal exponent a = 0.72 ±0.04, a prediction 
that should be experimentally testable. 

Whatever the microscopic origin of the novel FS state, 
general symmetry principles dictate that the long length 
scale phenomenology is described by an effective Landau- 
Ginzburg free energy functional 
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where ip is the superconducting order parameter and A 
is the electromagnetic vector potential with B = V x A. 
The constants m, a, b, K, t, and u are experimentally 
measurable phenomenological parameters, with a and t 
changing their signs at the superconducting and ferro- 
magnetic decoupled transition temperatures respectively, 
and 00 — hc/2e the elementary magnetic flux quantum. 

The above model displays rich phenomenology It 
has been shown in detail by Greenside, et al. that 
for a range of physically realistic parameters (e.g., Xl/S,, 
\Ll\/K/\t\ « 0(10), [| with the large Abrikosov ratio 
required for a robust mixed state, and large exchange K 
necessary to suppress the competing spiral phase, rela- 
tive to the SV state), systems described by Fgl, EqJ|, 
exhibit a thermodynamically stable phase consisting of 
a spontaneous (i.e. H = 0) vortex state, in which su- 
perconducting vortices are generated by the spontaneous 
magnetization of the ferromagnetic component. 

Deep within this state, both iP{t) and M(r) are large, 
and therefore the London approximation applies. We 
take M(r) = Mon(r) to be of approximately uniform 
magnitude Mq and fluctuating direction n(r). We also 
take |'0(r)| to be a constant ■00 7^ 0, everywhere except at 
the locations of the vortices (where -0 = 0) and allow its 
phase 0{y) to vary subject to the circulation condition 
that it increases by an amount 27r along a closed path 
enclosing a vortex The magnetic induction B(r) is 
confined to thin flux filaments defined by the zeros of ■0- 
Using the circulation condition and the minimizing above 
free energy Fql with respect to ipQ, 9 and A, we elim- 
inate V0 and B in favor of vortex line conformational 
degrees of freedom and obtain: 
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where we have dropped unimportant constant pieces. 
t(r) = J2n I '^^^JT^^i^ ~ ^nis)) is the directed vortex 
line density, with magnitude given by the local flux den- 
sity (in units of 0o) and with direction tangent to the 



local vorticity; r„(s) is the locus of the n-th vortex line, 
parameterized by an arc length s. The vortex interaction 
V{r) is approximately given by the standard expression 
V{r) = e-'-/^^/(47rAir), with A| = m0g/(167r3n^|V'o|^) 
Q, cut off at scales smaller than the coherence length ^. 

At low temperature T and in the clean limit the vortex 
lines freeze into a lattice (as has been observed Q), di- 
rected on average along a spontaneously chosen magneti- 
zation direction iio, which we choose to be the z axis, and 
with the lattice constant determined by the average mag- 
netization Mq. Vortex lattice fluctuations are described 
by the 2d phonon displacement vector field u{z, R„), de- 
fined by r„(z) = (Xn + Ua:{z,Rn),Yn + Wy(z,R„),z), 
where R„ = (Xn,Yn) spans the vortex lattice, and we 
chose s = z, and by a small 2d vector (5n _L z defined by 
n = (1 — (5n • i5n)^/^z -I- 6n. 

Integrating out Sn and going to the continuum, we ob- 
tain 1^,^ the effective elastic Hamiltonian: 
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where k is the vortex line curvature modulus and /i and A 
are the Lame coefficients, /i, A and k are all expressible 
in terms of the parameters of the original model Eqj|, 
and Uij = ^{diUj + djUi — dzUidzUj) is the rotationally 
invariant nonlinear strain tensor, with i,j G x,y. This 
form of Hei is strictly dictated by general symmetry con- 
siderations, specifically: the underlying rotational sym- 
metry of the SV solid about any arbitrary axis guarantees 
that only the fully rotationally invariant combination Uij , 
given above, appears and that the vortex line tilt modu- 
lus (C44) vanishes identically (but see below) 

Although the softness of the SV lattice suggests, via 
the Lindemann criterion, a lower melting temperature, 
our long wavelength description cannot directly predict 
the melting temperature, because thermal fiuctuations 
are dominated by short distance modes (near the Bril- 
louin zone boundary). Instead, here we focus on the 
much more interesting, robust and universal effects of 
quenched disorder, all of which (at long wavelengths) can 
be incorporated into our model by adding to Hei 
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The first term accounts for the random positional pin- 
ning (at reciprocal lattice vectors G) of the vortex lattice, 
with VG(r) a complex random pinning potential, whose 
statistics can be accurately represented by zero-mean, 
short-ranged, Gaussian correlations: Vg(i")V'q, (r') = 
Av(5QQ/(5'*(r — r') The second term reflects the pres- 
ence of an additional orientational, tilt disorder. The tilt 
held h(r) is a random vector with s/iort-ranged isotropic 
spatial correlations, whose statistics can also be taken 
zero-mean Gaussian, with K{Y)hj{v') — AS,,jS^{r - r') 
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0. In the case of magnetic impurities, it is the ran- 
dom local fields, acting on M, i.e., h(r) • n, which via 
the Higg's mechanism, i5n « dzU lead to the tilt disorder 
piece of Hdis- In the case of nonmagnetic impurities, it 
is the short-scale anisotropy in the spatial distribution of 
vortex line pinning centers, together with lack of rota- 
tional symmetry in the SV state, that leads to the tilt 
disorder in Eq.^, as can be easily illustrated by short- 
scale modes thinning procedure. ^ (As in conventional 
vortex lattices, in (i > 2, the pinning of the long-scale 
vortex density distortions is irrelevant at long scales). 

Adapting Larkin's Q calculation to the "soft" elastic- 
ity of the SV lattice we find that in d dimensions and 
at short length scales, mean-squared lattice distortions 
(in domains of size L±{^ L^)) grow as w^ms L^/^^"^. 
Hence, translational order is destroyed for d = 3 < 
d^^ = 9/2, implying (in the weak disorder limit) di- 
vergently anisotropic Larkin domains spanned by 3d 

Larkin lengths « {a^{KBiy/*/AvY^^ and « 

1/3 

(a^Ki?s/Ay) , with a the vortex spacing, s ^ L,T, and 
Bl = 2/x-l-A, Bt — Beyond these Larkin length scales, 
the Larkin approximation breaks down and the full non- 
linear nature of the random-field pinning potential must 
be taken into account, leading to much more weakly 
(specifically, logarithmically) divergent lattice fluctua- 
tions as in conventional vortex lattices |]l0|| . This weak 
divergence is completely overwhelmed for d < due = 7/2 
by the far stronger fluctuations induced by the tilt disor- 
der, which diverge algebraically: u^ms ~ (ignor- 
ing nonlinear elastic effects) out to arbitrarily long length 
scales 

Ignoring nonlinear elastic terms, we obtain the 
disorder-averaged phonon correlation function Cs{v) ~ 
Max[A/(K3s5)i/4r7/2-d^ A/(KS-'^Bf-i)i/2z7-2rf] X 

0(1). The spatial (Larkin-like) correlation lengths, 
measurable in neutron scattering experiments, are 
defined by C,(Min[a, C.]) 
U « [aH^'B^,)^/yAf^'^'-'\ i 
Hence, in 3d, in the presence of arbitrarily weak disorder 
A, on length scales longer than the translational 
order of SV solid decays rapidly, that is, it is short- 
ranged. This, however, does not imply that the system is 
equivalent to the SV liquid state, or to a fully disordered 
vortex glass as it would be if dislocation loops unbound. 
Using duality methods, introduced recently by two of 
the authors we have shown that a translationally 
disordered, but topologically ordered "columnar elastic 
glass" (CEG) phase is in fact possible in 3d. This is anal- 
ogous to our recent results in the context of randomly 
pinned smectic liquid crystals and columnar phase 

of discotic liquid crystals confined in aerogel, as well 
as the Bragg glass conjecture for conventional vortex 
lattices and the random field XY model 14|. 

In this topologically ordered phase the anharmonic 
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elastic terms in the Hamiltonian lead to drastically mod- 
ified elastic scaling on length scales longer than the 
nonlinear crossover length scales = Stt^k^/A and 
^NL — i.iNL)'^{Bs/ k)^^'^ ■ In renormaliztion group (RG) 
language, on scales longer than our system flows 

away from a zero-disorder Gaussian fixed point and ap- 
proaches a nontrivial T = anharmonic, disordered fixed 
point. This new fixed point can be studied using stan- 
dard RG analysis |^ , and we predict that at scales longer 
than out to an arbitrary length scale (but see be- 

low), the elastic moduli k, A, and the disorder variance 
A become singular functions of wavevector k |5| : 

«(k) = {k.^NLr'^ Uk±iNL/ikdNLf) , (6a) 
M(k) = MO [k.^NLT'' Mk±^NL/{k.eNL)'^) , (6b) 

A(k) = .9,Mo (fc^^L)"^ Mk±iNL/{kdNL)^) , (6c) 
A(k) = Ao (fc.e^j'-'"^ fAik^^hL/ik^^NL)'^) , (6d) 

where f^/^/Ai^) is independent of a; if a; ^ 1 and goes 
like X to the power —rj^/C,, rj^/C, and (1 — 2r7„)/C, re- 
spectively, when j; 3> 1, so that, e.g., K(k) cx k'^"^^ 
for fcj_^^^ 3> (fcz'^Afi)^- kq, Mo and Aq are the bare 
values of the elastic constants and the anomalous expo- 
nents are universal, and to leading order in e = \~ d are 
given by r?^ = 1.478e = 0.74 and r^^ = 0.6919e = 0.346, 
where the second equalities are their values in 3d. The 
constant is also universal, and to leading order in 
e g^, = —0.03272 -t- 0(e). The anisotropy exponent 
C = 2 — (77k-|-?7^)/2, and therefore the two independent ex- 
ponents and 77^ completely characterize the anomalous 
elasticity, with K(k) and A(k) diverging with linearly- 
related exponents, and /i(k) and A(k) both vanishing in 
exactly the same way as k — > 00. The negative universal 
amplitude ratio g^ implies that the SV glass exhibits a 
negative universal Poisson ratio Up = —0.0338 + 0(e). 

For length scales longer than ^tvl this strong, power- 
law anomalous elasticity will alter the behavior of all 
physical observables of the SV solid, such as, e.g., the 
width of the structure function peak S'(k) measured in 
neutron scattering and the behavior in its tails. It also 
leads to a non-Hookean, i.e., nonlinear, stress-strain re- 
lation even for arbitrarily small applied stress cr. 

To see this, consider a purely compressive stress, aij — 
aSij , which adds a term cr V • u = auu -\- ^IS^up to the 
Hamiltonian. In Fourier space, the second, symmetry 
breaking, term becomes crfc^ |u(k)p and begins to dom- 
inate over the K(k)fc^|u(k)p term once crfc^ > K(k)fc^. 
This clearly happens for kz's less than a critical kc 
given by crfc^ — K{kz = fcc,k_L = 0)/c^. Using Eq.|6a| 
for K(k) in this expression and solving for kc, we find 

kci<7) = (a/«o)^/('-''»'(?^L)''^/^'-'"'^- 

Now, for sufficiently weak stress cr, I/'Catl ^ kc{(j), the 
stress-induced rotational symmetry breaking term is sub- 
dominant to the vortex curvature energy K(k)fcf |u(k)p 
for wavevectors kz ^ kc{a). Consequently, on interme- 
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diate length scales, kc{<j) < < kf^i^, the SV glass ex- 
hibits anomalous elasticity summarized by Eqs.^. How- 
ever, on longer length scales, fcz ^ fcc(f), the apphed 
stress energy dominates, suppressing fluctuations and 
cutting ofl^ the anomalous elasticity at kc{<j). There- 
fore, the elastic moduli saturate^ for k^ <C kc{a), at 
their values at k^ — kda), fc^ = given by Eq.||, i.e., 
fi{k 0) ^ n{k, = kc,ki_ = 0) oc /c?' and A(k ^ 0) ^ 
X{kz ~ kc, kj^ = 0) cx kc^ . This implies that the Young's 
modulus Y{a) = 4^(^ + A)/(2^ + A) = 4/xo(l + g*)/(2 + 

9*){kceNL)''' « 2Aio(fccC^L)''- = 2^io{<j/'JNLf, with 

a < CTjvL, /3 = r/p/(2 - 77„), and (TArL = koCjvl- Hence, 
we find for stress a < cttvl, a nonlinear (non-Hookean) 
strain-stress relation: eda) = V • u = (T/y(cr) oc cr", 
with a = 1 - /3 = 1 - r/^/(2 - 77,) w 0.72. 

Now, it is easy to see that H = Hz, applied along the 
direction of the SV lattice, acts as a compressive stress, 
a = -ffS(O), and the induced SB{H) = B{H) - B{0) 
plays the role of the compressive strain, ec = SB/B{0). 
Hence H allows us to directly probe the anomalous elas- 
ticity, with the non-Hookean elasticity leading to the 
nonlinear and universal B{H) relation, for H < Hnl — 
crNL/B{0), given in Eq.|l| and Fig.Q. 

Of course, as with all crystalline ferromagnets, crys- 
talline symmetry breaking fields will explicitly break ro- 
tational invariance, even in the absence of H. They 
will lead to a non-zero tilt modulus, i.e., a term 
iVcr. |9zUp in the Hamiltonian, which will cut off the 
anomalous elasticity on length scales longer than = 

For weak Kr., our predic- 
tions will apply over a wide range of length scales Lz^± 
satisfying ^^^^ » Lz,_l ^ ^nl- On longer length scales 
the SV solid elasticity will cross over to the conventional 
"tension" elasticity of vortex lattices |^ , exhibiting a lin- 
ear relation SB cx H, for H < Her. = Vcr./B{0). 

The crystal anisotropy also leads to metastability for 
\H\ < Her., thereby allowing experimental studies of H 
applied in the direction opposite to that of M. In the 
limit \H\ ^ Her., the vortex solid becomes "soft" again, 
with the crystal-field generated tilt modulus vanishing as 
Her. ~ I^^Ij and the system exhibiting anomalous elastic- 
ity out to arbitrary length scale at this finely field-tuned 
point. However, because the rotational symmetry is nev- 
ertheless broken (with restoration only at the quadratic 
(dzu)^ level), we expect that the universality class and 
therefore the anomalous exponents will be distinct from 
that of the CEG studied here. || 

As argued in the case of vanishing crystal anisotropy, 
here too, the anomalous elasticity leads to a non-trivial 
field dependence of the Young's modulus. We find 
Y{H) oc {Her. ~ \H\)^ , with /3' > a universal exponent 
analogous to, but distinct from the f3 of the CEG. How- 
ever, in contrast to the rotationally invariant case, here 
the flux density SBiH) = H/Y{H) = H{Her - \H\)-^' 
diverges as approaches Her. from below. This signals 
a breakdown of the analysis, which assumes 6B{H) <C 



Bifi), and suggests that the metastability limit occurs at 
fields Hms determined by 5B{Hms) ~ B{Q). 

Finally, we note that ErNi2B2C has a strong uniax- 
ial anisotropy, with M preferentially lying in a nearly 
isotropic xz-plane. The appropriate model in this case 
resembles that given by Eq.^, except that in the strain 
tensor Uy , Ui Ux, and Uy has a conventional quadratic, 
tension elasticity. A similar model has recently been con- 
sidered by us in the context of the columnar phase of dis- 
cotics confined in an anisotropic random medium, such 
as, e.g., strained aerogel ||l2|. The resulting "m = 1 elas- 
tic glass" phase is also topologically stable with anoma- 
lous elasticity, and has additional interesting features, 
most notably, that it exhibits short-ranged translational 
order in one direction {x) and gwasi-long-ranged trans- 
lational order in the other (y), a prediction that should 
be testable via neutron scattering and decoration exper- 
iments. We conclude by noting that the "soft" elastic- 
ity of the SV solid should also manifest itself in many 
other physical observables, most notably in dynamics, 
exhibiting novel current-voltage characteristics and en- 
hanced critical currents, sensitive to and with nontrivial 
dependence on H. 
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